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Abstract 

All presently available results lead to the conclusion that nonextensivity, 
in the sense of nonextensive statistical mechanics (i.e., q ^ 1), does not modify 
anything to the second principle of thermodynamics, which therefore holds in 
the usual way. Moreover, some claims in the literature that this principle can 
be violated for specific anomalous systems (e.g., granular materials) can be 
shown to be fallacious. One recent such example is analyzed, and it is sug- 
gested how q ^ 1 distributions could in fact restore the validity of macroscopic 
time irreversibility, a cornerstone of our present understanding of nature. 

The Second Principle of Thermodynamics, i.e., time irreversibility of the macroscopic 
world, appears to be one of the most solid laws of theoretical physics. Usual statistical me- 
chanics, i.e., Boltzmann-Gibbs statistical mechanics, is consistent with this principle through 
the celebrated if -theorem. We address here what happens with this principle within the 
nonextensive statistical mechanics proposed in 1988 [|]] and characterized by an entropic 
index q (the particular case q = 1 recovers standard Boltzmann-Gibbs statistical mechan- 
ics). The answer seems to be very simple: no hope for moto perpetuo within this formalism! 
Indeed, the if-theorem appears to be g-invariant, more precisely, excepting for very quick 
microscopic fluctuations, S q cannot decrease (increase) with time if q > (q < 0), a fact 
which is consistent with the concavity (convexity) of S q , with regard to the set of proba- 
bilities, for q > (q < 0) (hence, S q is respectively maximal and minimal for q > and 



q < 0). This g-invariance was first pointed by Mariz || assuming detailed balance, and 
since then by many other authors (see [Q] and references therein) in a variety of generic 
situations. In addition to this, the present nonextensive thermostatistics provides a possible 
path for solving amazing "violations" of the second principle of thermodynamics. The rest 
of this paper is dedicated to the re-analysis of one such example that appeared recently in 
the literature of granular materials. 

Through molecular dynamics simulations using the inelastic hard sphere model to mimic 
fluidized granular media, Soto, Mareschal and Risso addressed f§] a very interesting question, 
namely the validity (or not) of Fourier's law and of the second principle of thermodynamics 
for dissipative macroscopic systems. They simulated a simple d = 2 system in the presence of 
a "vertical" (y axis) gravitational field, the top and bottom edges of a square box providing, 
to the arriving particles, a Maxwellian distribution of velocities at a fixed (dimensionless) 
temperature T = 1. The average density no and the dissipation coefficient qsMR (we use 
here the notation qsMR instead of the original notation q [|J in order to avoid confusion with 
the nonextensive entropic index q) were held fixed during a given simulation {qsuR — 
corresponds to the conservative limit, i.e., elastic collisions). They measured, as functions 
of the height y, the density n, the temperature T and the heat flux J. After long runs the 
system achieved a stationary state, like the one illustrated in their Fig. 1. The authors 
verified that T(y) exhibits a minimum at y — y*, whereas J vanishes at a value of y larger 
than y*. Consequently, the authors conclude that Fourier's law, which reads J = —kVT, 
k being the thermal conductivity, is violated. Moreover, for all values of y not exceedingly 
above y*, the heat appears to flow from the cold regions to the hot regions, therefore the 
authors consistently conclude that the second law of thermodynamics is violated as well. The 
way out from this paradoxal situation that Soto et al advocate is that an extra, n-dependent, 
term must be included in the expression of J. 

The purpose of the present remark is to argue that the proof provided by Soto et al would 
be very neat were it not a subtle weakness that we now address. The quantities n and J 
observed in the simulations are calculated essentially through their definitions. This is not 
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the case of T, which the authors set equal to the average kinetic energy per particle. This 
equality holds for the well known Maxwellian distribution, an assumption which is certainly 
false inside the box (see, for instance, flU). In order to illustrate the subtleties, let us use 
nonextensive statistical mechanics. Indeed, the use of such generalized statistics for the 
present system is quite natural since the stationary distribution observed by Taguchi and 
Takayasu in a vibrating-bed fluidized state precisely corresponds, as we shall see, to q = 3 
(similarly to what occurs in Levy and correlated anomalous diffusions f7|), whereas they 
observe a Maxwellian distribution (i.e., q = 1) in the so-called solid state. The possibility 
of applicability of this formalism to such materials has also been recently mentioned by 
Herrmann ||. 

For the Hamiltonian Ji = | J2u=\ the nonextensive formalism yields, for the velocity 



distribution probability, 



Pq({Vi}) 



l-(l-q)(H-U q )/T q 



l-q 



(1) 



where Z q = J[nf =1 d Vi J [l - (1 - q)(H - U q )/T q \ U q = /[nf^dv^ H/ f[Uf =1 dVj 
and T q = TZ x ~ q . This equilibrium distribution has been obtained by optimizing S q 



P q q 



1 ~~ / n^Li^Vj p q J /(q — 1) with appropriate constraints (Boltzmann constant ks has been 
taken equal to unity). It recovers the usual Maxwellian distribution for q = 1 and a power- 
law for q ^ 1. Also, we verify that 1/T — dS q /dlf q . Finally, we can rewrite Eq. (1) as 
follows: 



Pq 



l-(l-q)H/r 



l-q 



(2) 



where Z' q = J[llf =1 dv j \ [l - (1 - q)H/T'\ ^ and V = T q + (1 - q)U q . 

We verify that, for q — l,T q — T' — T, but such simplification disappears if q ^ 1, which 
might be the case for the system considered by Soto et al. In particular, if we identify p q q {p q 
being the N — 1 distribution (2)) with the distribution numerically obtained in || for the 
vibrating-bed fluidized state, we obtain q/(q — 1) = 3/2, hence q = 3. In other words, some 



further considerations are necessary before drawing any conclusion about the validity (or not) 
of Fourier's law and of the second principle of thermodynamics. A remarkable simplification 
does exist due to the fact that Yamano and others (see QTUJ and references therein) have 
recently shown [[J that each quadratic velocity term contributes to U q with T q /2. But, for 
the rest, the entire Fourier law would have to be rededuced, possibly in the q ^ 1 scenario. 
A variety of mathematical tools |ll|] (e.g., the g-generalization of Kubo's linear response 
theory) are already available in the literature for performing such generalization. It is clear 
nevertheless that this is not a trivial task. Only after such an analysis, the interesting 
computer simulations of Soto et al could led us to further conclusions. It could even happen 
that, in the present granular matter ffuidized state, the generalized J vanishes precisely 
when T q attains its minimum! 
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